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Appendix A 

<?xmi version="1.0" encoding = "UTF-8" ?> 

- <model name = "FieldKorosNoyesModer> 

§ <notes> 

<hl>Field-Koros-Noyes Model of BZ Reaction</hl> 

- <table border="0" cellspacing = "0" cellpadding = M 2 n > 

- <thead> 

- <tr> 

10 <th align="left" vaiign="middle" bgcolor="#eeeeee">Citation</th> 

</tr> 
</thead> 

- <tbody> 

- <tr> 
15 <td> 

R.J. Field and R.M.Noyes,J.Chem.Phys.60,1877 (1974) ; 

R.J. Field, E.Koros,R.M.Noyes,JACS 94,8649 (1972);R.J. Field, R.M.Noyes,Nature 
237,390 (1972) This implementation is taken manufactured by J.D. Murray, 
"Mathematical Biology" (1989) page 181. 

20 <a href= ,,n /> 
</td> 
</tr> 
</tbody> 
</table> 

25 <table border= M 0" cellspacing = "0" cellpadding = "2"> 

- <thead> 

- <tr> 

<th align="left M valign = "middle" bgcoIor="#eeeeee">Description</th> 

</tr> 
30 </thead> 

- <tbody> 

- <tr> 

<td> Field Noyes Version of Belousov- Zhabotinsky Reaction. Br03 is held 

constant; HOBr is typically ignored, and can be replaced by an empty- set. The 
35 stoichiometry f is typically taken as 1/ 2 or 1 (denominator 1 or 2 in SBML) 

.</td> 
</tr> 
</tbody> 
</table> 

40 <table border="O n cellspacing = "0 M cellpadding = "2 M > 

- <thead> 

- <tr> 

<th align =" left" valign = ,, middle" bgcolor="#eeeeee">Rate constant </th> 
<th align="left" valign="middle" bgcolor="#eeeeee u >Reaction</th> 

45 </tr> 

</thead> 

- <tbody> 

- <tr> 

<td>kl = 1.3</td> 
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<td>Br + Br03 -> HBr02 + HOBr</td> 

</tr> 

- <tr> 

<td>k2 = 2000000</td> 
5 <td>Br + HBr02 -> HOBr^2</td> 

</tr> 

- <tr> 

<td>k3 = 34</td> 

<td>Br03 + HBr02 -> Ce^2 + HBr02^2</td> 

10 </tr> 

- <tr> 

<td>k4 = 3000</td> 

<td>HBr02 A 2 -> Br03 + HOBr</td> 

</tr> 
IS <tr> 

<td>k5 = 0.02</td> 
<td>Ce -> Br A f</td> 

</tr> 
</tbody> 
20 </table> 

- <table border="0" cellspacing= M 0" cellpadding="2 M > 

- <thead> 

- <tr> 

<th align="left" valign = M middle ,, bgcolor= M #eeeeee">Variable</th> 
25 <th align="left" valign="middle n bgcolor="#eeeeee">IC </th> 
<th align = "left" valign = n middle M bgcolor= n #eeeeee n >ODE</th> 

</tr> 

</thead> 

- <tbody> 
30 <tr> 

<td>Br</td> 
<td>0.003</td> 

<td>Br'[t] == -(kl*Br[t]*Br03[t]) + f*k5*Ce[t] - k2*Br[t]*HBr02[t]</td> 

</tr> 
35 <tr> 

<td>Ce</td> 
<td>0.05</td> 

<td>Ce"[t] == -(k5*Ce[t]) + 2*k3*Br03[t]*HBr02[t]</td> 

</tr> 
40 <tr> 

<td>HBr02</td> 
<td>0.001</td> 

<td>HBr02'[t] == kl*Br[t]*Br03[t] - k2*Br[t]*HBr02[t] + 
k3*Br03[t]*HBr02[t] - k4*HBr02[t] A 2</td> 
45 </tr> 

- <tr> 
<td>HOBr</td> 
<td>0</td> 

<td>HOBr'[t] == kl*Br[t]*Br03[t] + 2*k2*Br[t]*HBr02[t] 
50 + k4*HBr02[t] A 2</td> 

</tr> 
</tbody> 
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</table> 
</body> 
</notes> 

- <listOfCompartments> 

5 <compartment name="BZ" /> 
</listOfCompartments> 

- <MstOfSpecies> 

<specie name="Br" initialAmount="0-003" compartment="BZ" boundaryCondition = "false" 
/> 

10 <specie name= M Br03" initialAmount= M 0.1" compartment="BZ M boundaryCondition = "true n 
/> 

<specie name="Ce" initialAmount=' , 0.05 M compartment= n BZ M boundaryCondition= M false" 

/> 

<specie name="HBr02" initialAmount="0.001" compartment= n BZ" 
15 boundaryCondition="false M /> 

<specie name="HOBr" initialAmount="0 ,, compartment="BZ" boundaryCondition = "false" 
/> 

</HstOfSpecies> 

- <listOfReactions> 

29 <reaction name="Reactionl" reversible="false"> 

- <listOfReactants> 
<specieReference specie= ,, Br" /> 
<specieReference specie= , 'Br03" /> 

</listOfReactants> 
2§ <listOfProducts> 

<specieReference specie= ,, HBr02" /> 
<specieReference specie= ,, HOBr ,, /> 

</MstOfProducts> 

- <kineticLaw formuIa = "Br*Br03*kl"> 

30 <NstOfParameters> 

<parameter name= ,l kl" value="1.3" /> 

< / 1 i stOf Pa ra m ete rs > 
</kineticLaw> 
</reaction> 

3§ <reaction name= ,, Reaction2 ,, reversible="false"> 

- <MstOfReactants> 
<specieReference specie= ,, Br n /> 
<specieReference specie="HBr02 M /> 

</listOfReactants> 
40 <listOfProducts> 

<specieReference specie="HOBr M stoichiometry="2 M /> 
</listOfProducts> 

- <kineticLaw formula = ,, Br*HBr02*k2"> 

- <listOfParameters> 

45 <parameter name = "k2" value= 2000000 /> 

< / 1 i stOf Pa ra m ete rs > 
</kineticl_aw> 
</reaction> 

- <reaction name= n Reaction3" reversible="false"> 
50 <listOfReactants> 

<specieReference specie= ,, Br03 M /> 
<specieReference specie="HBr02' , /> 
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</listOfReactants> 

- <listOfProducts> 

<specieReference specie="Ce n stoichiometry="2" /> 
<specieReference specie="HBr02" stoichiometry="2 M /> 
5 </listOfProducts> 

- <kineticLaw formula = ,, Br03*HBr02*k3 M > 

- <listOfParameters> 

<parameter name="k3" value="34" /> 

< / 1 i st Of Pa ra m ete rs > 
10 </kineticLaw> 

</reaction> 

- <reaction name="Reaction4 M reversible="false M > 

- <NstOfReactants> 
<specieReference specie= M HBr02 n /> 

15 </IistOfReactants> 

- <listOfProducts> 
<specieReference specie="Br03 M /> 
<specieReference specie= ,, HOBr H /> 

</MstOfProducts> 
20 <kineticLaw formula = "HBr02 /v 2*k4"> 

- <NstOfParameters> 

<parameter name="k4" value="3000 n /> 

< / 1 i st Of Pa ra m e te rs > 
</kineticLaw> 

25 </reaction> 

- <reaction name="Reaction5" reversible="false"> 

- <listOfReactants> 
<specieReference specie="Ce n /> 

</MstOfReactants> 
30 <listOfProducts> 

<specieReference specie="Br" stoichiometry^'l" denominator= n 2" /> 
</NstOfProducts> 

- <kineticLaw formula = "Ce 5 * c k5 n > 

- <listOfParameters> 

35 <parameter name= n k5" value="0-02" /> 
<parameter name="f" value= ,, 0-5 n /> 

</NstOfParameters> 

</kineticLaw> 

</reaction> 
40 </NstOfReactions> 

</model> 

</sbmI> 



